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ABSTRACT 

We present deep, accurate and homogeneous multiband optical (U,B,V,I) photometry of the Carina 
dwarf spheroidal galaxy, based on more than 4,000 individual CCD images from three different ground- 
based telescopes. Special attention was given to the photometric calibration, and the precision for the 
B, V, and /bands is generally better than 0.01 mag. We have performed detailed comparisons in the 
V, B-V and V, B-I color-magnitude diagrams (CMDs) between Carina and three old, metal-poor 
Galactic Globular Clusters (GGCs, M53, M55, M79). We find that only the more metal-poor GCs 
(M55, [Fe/H]=-1.85; M53, [Fe/H]=-2.02 dex) provide a good match with the Carina giant branch. 
We have performed a similar comparison in the V, V-I CMD with three intermediate-age clusters 
(IACs) of the Small Magellanic Cloud (Kron 3, NGC 339, Lindsay 38). We find that the color extent 
of the SGB of the two more metal-rich IACs (Kron 3, [Fe/H]=-1.08; NGC339, [Fe/H]=-1.36 dex) 
is smaller than the range among Carina's intermediate-age stars. Moreover, the slope of the RGB 
of these two IACs is shallower than the slope of the Carina RGB. However, the ridge line of the 
more metal-poor IAC (Lindsay 38, [Fc/H]= 1.59 dex) agrees quite well with the Carina intermediate- 
age stars. These findings indicate that Carina's old stellar population is metal-poor and appears to 
have a limited spread in metallicity (A[Fe/H]=0.2-0.3 dex). The Carina's intermediate-age stellar 
population can hardly be more metal-rich than Lindsay 38 and its spread in metallicity also appears 
modest. We also find that the synthetic CMD constructed assuming a metallicity spread of 0.5 dex for 
the intermediate-age stellar component predicts evolutionary features not supported by observations. 
In particular, red clump stars should attain colors that are redder than red giant stars, but this is not 
seen. The above results are at odds with recent spectroscopic investigations suggesting that Carina 
stars cover a broad range in metallicity (A[Fe/H]~l-2 dex). We also present a new method to estimate 
the metallicity of complex stellar systems using the difference in color between the red clump and the 
middle of the RR Lyrae instability strip. The observed colors of Carina's evolved stars indicate a 
metallicity of [Fe/H]=-1.70±0.19 dex, which agrees quite well with spectroscopic measurements. 
Subject headings: galaxies: dwarf — galaxies: stellar content — stars: evolution — stars: fundamental 
parameters 
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1. INTRODUCTION 

Dwarf galaxies play a fundamental role in several as- 
trophysical problems. Current cosmological simulations 
predict dwarf satellite populations significantly larger 
than the number of dwarfs observed near giant spirals 
like the Milky Way and M31. This discrepancy is called 
the "missing satellite problem" and is a challenge to the 
currently most popular cosmological model: the A Cold 
Dark Matter paradigm (Klypin et al. 1999; Moore et al. 
2006; Madau et al. 2008). However, it has not yet been 
established whether this discrepancy is due to limitations 
in the theoretical modeling or to observational bias at the 
faint end of the dwarf galaxy luminosity function (Kli- 
mentowski et al. 2009; Koposov et al. 2009; Kravtsov 
2010). 

Moreover, detailed analyses of the distribution of dwarf 
galaxies in multi-dimensional parameter space indicate 
that they follow very tight relations. In particular, Prada 
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& Burkert (2002) suggested that Local Group (LG) dwarf 
galaxies show strong correlations (Fundamental Line - 
FL) among mass-to- light (M/L) ratio, surface brightness, 
and metallicity. They explained the correlation between 
M/L ratio and metallicity using a simple chemical enrich- 
ment model in which the hot metal-enriched gas, at the 
end of the star-formation epoch, is lost via continuous 
galactic winds (Mac Low & Ferrara 1999). In a recent 
investigation Woo & Chiueh (2009) reached similar con- 
clusion, i.e., LG dwarf galaxies follow a one-parameter re- 
lation driven by the total stellar mass. They also found 
that dwarf spheroidals (dSphs) appear to be, at fixed 
stellar mass, systematically more metal-rich than dwarf 
irregulars (dls, see also Mateo (1998)). Futhermore, the 
FL of dls in a five-dimensional parameter (total space 
mass, surface brightness, rotation velocity, metallicity, 
star formation rate) is linear and tight, and extends the 
scaling relations of giant late-type galaxies into the low- 
mass regime. On the other hand, the FL of dSphs as de- 
fined in a four-dimensional parameter space (total mass, 
surface brightness, rotation velocity, metallicity) is also 
linear and very tight, but it does not lie on an extrapo- 
lation of the scaling relations of giant early-type galax- 
ies. Iron and heavy-element abundances, the spread in 
these quantities, and the star-formation history (SFH) of 
LG dwarf galaxies, are crucial observational constraints 
on the physical mechanisms responsible for the scaling 
relations described above, and on their implications for 
different cosmological models (Orban et al. 2008). Spec- 
troscopic investigations based on high-resolution spectra 
of bright red giants (RGs) in several LG dSphs indi- 
cate a large spread in iron abundance (Shetrone et al. 
2001, 2003). This is also supported by metallicity es- 
timates from the Call triplet method (Battaglia et al. 
2008). In this context, the Carina dSph is particularly 
relevant: (i) It is relatively close and its central density is 
modest (p=0.17M Q /pc 3 , Mateo (1998)) so it is easily re- 
solved with ground-based telescopes (Mould & Aaronson 
1983). (ii) It is the LG dSph most clearly showing mul- 
tiple, widely separated star formation episodes (Mighell 
(1990, 1997); Smecker-Hane et al. (1994, 1996); Hurley- 
Keller et al. (1998); Hernandez et al. (2000); Harbeck et 
al. (2001); Dolphin (2002); Rizzi et al. (2003)) with a 
tail that might extend to less than one Gyr ago (Monclli 
et al. (2003)). (Hi) It is also a benchmark to constrain 
the pulsation properties of old (RR Lyrae, Saha et al. 
(1986)) and intermediate-age (dwarf Ccpheids, Mateo et 
al. (1998); Anomalous Cepheids, DalPOra et al. (2003)) 
variable stars, (iv) High-resolution spectra are available 
for a sample of ten bright RGs. The mean metallicity 
is [Fe/H]=-1.69 with cr=0.51 dex (Koch et al. 2008). 
Independent measurements by Shetrone et al. (2003), 
using high resolution spectra for five bright RGs, pro- 
vided a mean metallicity of [Fe/H]=-1.64 and cr=0.2 dex. 
(iv) Medium-resolution calcium triplet measurements are 
also available for a large sample (437) of RG stars (Koch 
et al. 2006); their metallicity distribution shows a peak 
at [Fe/H]=-1.72±0.01 (metallicity scale by Carretta & 
Gratton 1997) and the metallicity ranges from ~-2.5 to 
~-0.5 dex. Using the same spectra, but different se- 
lection criteria (364 candidate Carina stars), Helmi et 
al. (2006) found the same peak in the metallicity distri- 
bution ([Fe/H]=-1.7±0.1 dex) and metallicities ranging 
from ~-2.3 to ~-1.3 dex (see their Fig. 2). 



Deep, accurate photometric investigations indicate an 
old (13 Gyr) stellar population and several distinct 
intermediate-age populations (2-6 Gyr) together with a 
blue plume of younger stars. The occurrence of both old 
and intermediate-age populations is also indicated by two 
distinct samples of core helium burning stars: an old, 
wide HB and a separate red clump (RC). The shape and 
thickness of the RGB (AB-V~0.02 for 20.5<U<20.75 
mag) suggest that the old and intermediate-age stellar 
populations have a minimal spread in chemical compo- 
sition. This seems inconsistent with the spectroscopic 
observations that indicate a metallicity distribution pos- 
sibly reaching extreme values near -3.1 and +0.1 dex 
(Zinn & West 1984, scale), or near -2.8 and -0.2 dex 
(Carretta & Gratton 1997, scale). 

Despite the criticisms raised by Koch et al. (2006), 
conclusions based on purely photometric indicators are 
hampered by only two major factors: (i) A decrease in 
photometric precision when moving from bright to faint 
RG stars might produce a spread in color mimicking a 
spread in metal content. In particular, Carina's stellar 
density is very low, requiring that photometry be ob- 
tained with multi-chip CCD cameras or with many dif- 
ferent pointings of a single CCD. A high degree of con- 
sistency in the photometric zero-points across the field 
is essential to avoid spurious broadening of the RG and 
Main Sequence (MS) loci, (ii) Due to its low surface 
density, large extent, and low Galactic latitude (-22°), 
the Carina sample is contaminated by foreground field 
stars. Their colors, unfortunately, are similar to the Ca- 
rina RGB, once again mimicking a broadening in color 
due to a spread in age/metallicity. 

Our group is involved on a long-term investigation of 
the stellar populations in the Carina dSph. In particu- 
lar, we have assembled our own and archival UBVI im- 
ages covering the entire body of the galaxy. In this in- 
vestigation we focus on comparing Carina with old and 
intermediate-age template clusters. We also present a 
new method for estimating the metal content of complex 
stellar systems from the difference in color between the 
middle of the RR Lyrae instability strip and the peak of 
RC stars. 

2. OBSERVATIONS AND DATA REDUCTION 

The data were collected in various observing runs be- 
tween 1992 December and 2005 January. They include 
images from three telescopes: the CTIO 1.5m telescope 
with a single Tektronix2K CCD, the CTIO 4m Blanco 
telescope with the MosaicII camera, and the ESO/MPG 
2.2m telescope with the Wide Field Imager camera (both 
proprietary and archival data). Details of the observa- 
tions will be presented in a future paper. The data pre- 
sented here represent 4,152 individual CCD images with 
essentially complete coverage of the central regions of Ca- 
rina (« 40' x 55', lacking only those regions obliterated 
by bright foreground stars), and some sampling of the 
galaxy's halo out to an extreme radius of ^108'. They 
were obtained in four photometric bands (B and / with 
the 1.5m telescope, and UBVI with both the 4m and 
2.2m telescopes). 

The data were reduced using the 
DAOPHOT / ALLFRAME package (Stetson 1987, 
1994). Individual PSFs were produced for each chip 
of every exposure, using semi-automated routines to 
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select bright, isolated PSF stars. Subsets of the data 
from each observing run were first reduced separately. 
Finally, everything was merged together and a single run 
of ALLFRAME was adopted to reduce all the images 
of the center of Carina, with separate reductions for 
non-overlapping outlying fields. Because of the multiple 
chips and pointings required to cover the galaxy, no 
star appears in every image; any given star may have 
up to 17 calibrated measurements in U, 156 in B, 
207 in V, and 70 in I. A total of 205,338 individual 
stars were catalogued and measured. Among these, 
72,595 have photometric measurements in all four filters, 
and 129,230 have at least V and either B-V or V-I. 
The remainder, cither extremely faint or located near 
the periphery of our coverage, have astrometry and 
instrumental magnitudes only. 

The Carina data were contained within 206 individual 
datasets, where a dataset is essentially the totality of 
data obtained with one CCD from one night of observing. 
These, along with 1,331 datasets from other nights and 
other telescopes, were calibrated to the current version 
of Stetson's (2000) photometric standard system, which 
is believed to be equivalent to that of (Landolt 1992) to 
well under 0.01 mag in each of B, V, and / (see also 
Stetson (2005)). The U photometric bandpass is more 
problematic, and in this paper we will employ the U- 
band data only for qualitative and relative comparisons. 

3. RESULTS AND DISCUSSIONS 

For a robust discrimination of candidate Carina mem- 
bers from foreground field stars and background galaxies 
we adopted the U-V,B-I color-color diagram (C-CD). 
The top panel of Fig. 1 shows that this is effective for 
distinguishing candidate Carina and field stars. Carina 
RGs (U-V>0.6, B-I>\A) attain, at fixed B-I, bluer 
U-V colors than field stars. The difference is mainly 
caused by a difference in the mean metallicity between 
Carina and the more metal-rich Galactic disk stars, but 
the difference in surface gravity between galaxy giants 
and foreground dwarfs may also play a role. In con- 
trast, hot horizontal-branch (HB) stars and relatively 
young MS stars in Carina (B-I<0A) are relatively free 
from field star contamination, since they are bluer than 
the thick disk and halo turnoffs; compact background 
galaxies represent the principal source of contamination 
among these bluer objects. Fortunately, this C-CD is 
also a good diagnostic to separate candidate Carina stars 
from blue background galaxies since the latter show, at 
fixed B-I, systematically bluer U-V colors than stars 
of any luminosity class or metallicity (apart from white 
dwarfs, which are not an issue here). 

Based on the above empirical evidence, we performed 
a series of tests to identify candidate Carina stars us- 
ing various color-magnitude diagrams (CMDs) together 
with the C-CD. The solid-outlined box in the top panel 
of Fig. 1 shows the final selection, while the bottom pan- 
els show the V,B-V CMD of the total sample (left), 
the candidate Carina stars (middle), and the candidate 
non- member objects (right). Note that after the selec- 
tion there remain some field stars with colors similar to 
Carina's RGB ( V<21.5, B-F~0.4-0.6; middle panel). 
However, RG stars are distributed along a very narrow 
sequence over the entire magnitude range (18<y<23), 
and the number of field stars lying near that RGB in both 



color and magnitude is appreciably reduced. The same 
statement applies to the intermediate-age red clump 
(RC: y-20.5, B-V~0M) and to the old HB (y-20.75, 
B-V ^0.65). These evolutionary phases are character- 
ized by very narrow distributions in either magnitude 
or color or both. The improvement compared to data 
already available in the literature is due to smaller pho- 
tometric uncertainties (ag_y<0.02 for V<21, see error 
bars in Fig. 1). Plain evolutionary arguments (Monclli 
ct al. 2003; Zoccali et al. 2003) suggest that small disper- 
sions in magnitude and color in these evolutionary phases 
imply a negligible spread in chemical composition. This 
finding supports the results based on B,V,I photometry 
of Carina RG stars obtained by Rizzi et al. (2003). The 
presence of a few stars in the right panel with magnitudes 
and colors typical of red HB (V~20.75, B-V~0A) and 
faint subgiant branch (SGB) stars (V~23, B-V~0.5) 
demonstrates that our C-CD selection is imperfect — 
especially at intermediate colors — and some real Carina 
members have probably been erroneously rejected. How- 
ever, in these photometric selections it is more important 
to keep probable non-members out than to keep possible 
members in (Calamida et al. 2009). 

To provide robust constraints on Carina's stellar con- 
tent, we compared it with stellar systems character- 
ized by different ages and metal abundances. In par- 
ticular, to constrain the old stellar population we se- 
lected three globular clusters (GCs)— M79 = NGC 1904, 
M55 = NGC 6809, and M53 = NGC 5024— with metal 
abundances ranging from -1.64±0.15 to -2.02±0.15, and 
low foreground reddenings (E(B-V)<0.09, see Table 1). 
The optical BVI photometry of these GCs is also on 
the Stetson (2005) system. The accuracy of the pho- 
tometry is certainly better than our knowledge of the 
foreground reddening toward these clusters and Carina. 
The V, B-V (left panels) and V, B-I (middle pan- 
els) CMDs of the selected GCs are shown in Fig. 2, to- 
gether with their individual ridge lines (red curves). We 
also selected three intermediate-age clusters (IACs) be- 
longing to the Small Magellanic Cloud (SMC)— Kron 3, 
NGC 339, Lindsay 38 — with similar ages and iron con- 
tents ranging from -1.08±0.12 to -1.59±0.10 dex (Glatt 
ct al. 2009, see Table 1). The photometry for these clus- 
ters is from ACS@HST data transformed into the V, I 
Johnson-Cousins bands using relations by Siranni et al. 
(2005). These clusters are not included in the homoge- 
nous photometry project (the original sources of the data 
are listed in Table 1), so we are unable to independently 
confirm the precision of the photometry or the accuracy 
of the calibrations. Data plotted in the right panels of 
Fig. 2 display the V, V-I CMDs for these clusters to- 
gether with their ridge lines and the outlines of their 
RCs. 

To compare the GCs with the old stellar populations in 
Carina, we arbitrarially shifted their ridge lines in mag- 
nitude and color to provide a good fit with the SGB 
(23.25<y<23.75, 0.45<5-V<0.60) in Carina. Then we 
inferred the difference in distance and in reddening be- 
tween Carina (/z=20.15, E(5-V)=0.03; Dall'Ora et al. 
(2003); Pictrzynski et al. (2009)) and the individual GCs. 
The distances and the cluster reddenings we found fol- 
lowing this approach agree, within the errors, with sim- 
ilar estimates available in the literature (see columns 4 
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and 5 in Table 1). The top and the middle left panels of 
Fig. 3 show that the GCs M79 and M55 provide a poor 
match to Carina, since hot HB stars are systematically 
brighter than Carina HB stars. On the other hand, both 
the GB and the HB of the more metal-poor GC (M53) 
match Carina quite well. The same outcomes apply if 
the comparison is performed using the V, B-I CMDs 
(see the middle panels in Fig. 3). The extinction in the 
BVI bands was estimated using the empirical relations 
provided by McCall (2004). 

The same approach was followed to compare the IACs 
with Carina stellar populations. The ridge lines were 
shifted in magnitude and color to provide a good fit to the 
lower envelope of intermediate-age SGB (22.5<V r <23.0, 
0.45< F-/<0.80) stars in Carina. The distances that 
we estimate are larger than those estimated by Glatt 
et al. (2009) (see Table 1) using isochrones. However, 
the current distances agree, within the errors, with SMC 
distances available in the literature (Romaniello et al. 
2009). The same outcome applies for the cluster redden- 
ings (see the right panels in Fig. 3). This comparison 
between Carina and the IACs revealed the following: (i) 
the color extents of the SGBs of the the two more metal- 
rich clusters (Kron 3 and NGC 339; i.e., the length of the 
horizontal sequence between the upper MS and the lower 
GB) are smaller than the color range covered by Carina's 
intermediate-age SGB stars; (ii) the slope of the RGB in 
the two more metal-rich clusters (Kron 3, NGC 339) is 
shallower than the observed slope of Carina's RGB; (Hi) 
the ridge line of Lindsay 38 is a fairly good match to Ca- 
rina's intermediate- age population, (iv) the RC contour 
of Kron 3 agrees with Carina's RC, while for the other 
two clusters the clumps are either slightly redder and 
brighter (NGC 339) or marginally fainter (Lindsay 38). 

These comparisons indicate that Carina's old popula- 
tion is relatively metal-poor, and its spread in metallicity 
is modest: at most 0.2-0.3 dex. Carina's intermediate- 
age population seems to have the same metallicity as 
Lindsay 38, or perhaps slightly lower. The spread in 
metal abundance in this stellar component also seems 
very limited, since the RGBs of the two more metal- 
rich IACs are considerably redder than the red envelope 
of Carina's RGB. The lack of accurate photometry for 
metal-poor ([Fe/H]<-1.6 dex) IACs with well populated 
RGBs prevents us from providing more firm constraints 
on Carina's intermediate-age population. These findings 
do not depend on the adopted metallicity scale, and if 
we adopt instead iron abundances on the Carretta et al. 
(2009) metallicity scale (see column 3 in Table 1), the 
metallicity range inferred for Carina is very similar. 

For a new estimate of the metal abundance of Carina 
we decided to use a new indicator. Evolutionary pre- 
scriptions for old and intermediate-age stellar structures 
indicate that the difference in color between the center 
of the RR Lyrae instability strip and the RC stars is 
strongly correlated with iron abundance. The correla- 
tion is caused by the mild dependence of the color of the 
instability strip on iron abundance and by the significant 
dependence of the color of RC stars when moving from 
metal-poor to metal- rich compositions (see Fig. 4). We 
estimated the predicted difference in color between RC 
stars (at the beginning of central helium burning) and the 
center of RR Lyrae strip (logX' e =3.85) from the large set 
of scaled-solar evolutionary models (Pietrinferni et al. 



2004) available in the BaSTI database 19 . We adopted 
evolutionary prescriptions that cover a wide range of 
iron abundances (-2.27<[Fe/H]<+0.06 dex), and an age 
range of 1-6 Gyr for the intermediate-age population 
and 13 Gyr for the old population. Note that the in- 
stability strip is minimally affected by cluster age when 
moving from 7-8 to 12-13 Gyr. In this context two 
robust evolutionary predictions need to be underlined. 
(i) When moving from metal-poor to metal-rich compo- 
sitions, the intermediate-age RC becomes redder than 
the old RGB for [Fe/H]=-1.79 dex (see the top panel 
of Fig. 4). Current evolutionary models predict that a 
change in metallicity from [Fe/H]=-1.96 to [Fe/H]=-1.49 
produces a change in V magnitude of RC stars of 0.03 
mag, but a change of 0.12 mag in the B-V color, (ii) 
The spread in magnitude of HB stars at the mean color 
of the RR Lyrae strip increases, as expected, with the 
spread in metallicity. Note that a change in metallic- 
ity from [Fe/H]=-1.96 to [Fe/H]=-1.49 implies a change 
in V magnitude of HB stars of 0.10 mag, but a change 
smaller than 0.01 mag in the B-V color of the instability 
strip (see the bottom panel of Fig. 4). 

Finally, we performed a linear regression between 
metallicity and the RC-HB difference in three different 
colors: 

[Fe/H] = (-3.02 ±0.03) + (3.89 ± 0.05)A(fl-V)j& 
[Fe/H] = (-3.47 ±0.04) + (2.37 ± 0.04) A(B-/)|g 
[Fe/H] = (-4.13 ±0.08) + (6.02 ±0.15)A(V-7)§^ 
The above relations (see also right panels in Fig. 5) 
indicate that the difference in the B-I color is signif- 
icantly more sensitive to [Fc/H] than the B-V color. 
The V-I is the least sensitive of all and shows a stronger 
nonlinear trend in the metal-poor ([Fe/H]<— 1.5) regime. 
We estimated the difference in color between the RC 
and the middle of the RR Lyrae instability strip (see the 
left panels in Fig. 5). We found AS- Vg<£=0.37±0.05, 
A5-/fg=0.75±0.05 and A V-J$£=0.40±0.05 mag, 
where the errors account for uncertainties in the pho- 
tometry and in the estimate of the mean color of the 
RR Lyrae strip. We applied the above relations and 
we found the following mean metallicities: [Fe/H]=- 
1.60±0.24(Afl-Vg%), [Fe/H]=-1.70±0.19(A5-/£<£) 
and [Fe/H]=-1.76±0.44(Ay-/|^). We adopted the 
estimate based on B-I, since it is the most precise. 

As a further test to constrain the spread in metallic- 
ity of the Carina sub-populations, we computed a se- 
ries of synthetic CMDs changing the spread both in age 
and in metal content. The set of evolutionary models 
used for these numerical simulations are the same used 
to derive the new metallicity indicator (Pietrinferni et 
al. 2004). The code adopted to compute the synthetic 
CMDs has already been discussed by Pietrinferni et al. 
(2004) and by Cordier et al. (2007). We only briefly 
mention the initial parameters adopted to compute the 
synthetic CMDs. The Initial Mass Function was mod- 
eled using a power law with a Salpeter exponent (a=- 
2.35), and we assumed a fraction of unresolved binaries 
of the order of 10% with a minimum mass ratio for the 
binary systems equal to 0.7. The synthetic CMDs were 
constructed assuming two star formation events, with a 

19 Evolutionary models can be download from http://www.oa- 
teramo.inaf.it/BASTI 
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mean age of 2.5-4 Gyr for the intermediate-mass and of 
11-13 Gyr for the low-mass sub-population (Monelli et 
al. 2003). We also assumed that the two star formation 
episodes include the same fraction of stars. Moreover, 
we assumed mean metallicities of [Fe/H]=-1.5 and of - 
2.0 dex for the intermediate-age and the old burst, re- 
spectively. In accordance with the above findings, we 
assumed the same spread in metallicity — ±0.1 dex — for 
the intermediate and the old star formation episode. Fi- 
nally, to mimic actual observations, synthetic photomet- 
ric errors were added using a Gaussian in the B and V 
bands with a dispersion equal to 0.02 mag. The above 
assumptions adopted to simulate the Carina observations 
are crude, but a detailed analysis is beyond the aims of 
this investigation. We plan to study the Carina star for- 
mation history in a forthcoming paper (Monelli et al. 
2010). Data plotted in the top panel of Fig. 6 show 
that the aforementioned assumptions provide evolution- 
ary features that are in reasonable agreement with the 
observations (see Fig. 2). In particular, both RC and 
HB stars are clearly separated in magnitude and in color 
from each other and from RG stars. Moreover, we se- 
lected two bins along the RGB for 20.75<m t ,<21.25 and 
21.50<m t ,<22.00 and we found that the spread in B-V 
color is ~0.06 mag. This value is quite similar, within 
the errors, to the spread in B-V color of the observed 
CMD at the same magnitude intervals, i.e., 0.06 mag. 
To constrain the impact of a spread in metal content, 
we used the same evolutionary ingredients adopted to 
construct the above synthetic CMD, but we assumed for 
each of the two stellar components a spread in metallicity 
of 0.2 dex. The spread in B-V color at the same mag- 
nitude intervals is slightly larger than observed, namely 
0.08 mag. 

Finally, we adopted for the intermediate-age stellar 
component alone a spread in metallicity of 0.5 dex. Data 
plotted in the bottom panel of Fig. 6 clearly show that 



such a spread in metallicity causes the RC stars to attain 
colors that are redder than the RG stars. The spread in 
B-V color at the same magnitude intervals is a factor of 
two larger than observed, i.e., 0.12 mag. Moreover, the 
synthetic CMD suggests the occurrence of subgiant and 
giant branch stars with colors that are redder than the 
subgiant of the old sub-population. The observed CMD 
docs not show these features. 

The comparison with empirical calibrators further sup- 
ports the contention that Carina hosts predominantly 
metal-poor stellar populations, in agreement with the 
spectroscopic investigations. However, the limited range 
in color covered by the RC stars appears incompati- 
ble with the spectroscopic results: the spread in iron 
abundance of the intermediate-age population is either 
small, or is counterbalanced by some other — as yet 
unrecognized — variable. The same conclusion applies 
for the old population, based on the limited magnitude 
spread among the old HB stars. These conclusions are 
also supported by synthetic CMDs constructed assuming 
a different spread in metallicity between the intermediate 
and the old sub-populations. The simulations indicate a 
spread in metallicity of the two sub-populations smaller 
than 0.2 dex. 
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TABLE 1 

Intrinsic parameters of old and intermediate-age calibrating clusters. 



Name 




[Fe/H] 


[Fe/H] c 




E(B-V) 


Age 


NGC 1904 - 


M79 


-1.64±0.10 a 


-1.58±0.02 


15.61±0.10 d 


0.01±0.02 d 


12.6±1.3 f 


NGC 6809 - 


M55 


-1.85±0.10 a 


-1.93±0.02 


13.69±0.10 d 


0.09±0.02 d 


12.4±1.7 f 


NGC 5024 - 


M53 


-2.02±0.15 a 


-2.06±0.09 


16.24±0.10 d 


0.01±0.02 d 


11.2±1.8 f 


Kron 3 




-1.08±0.12 b 


-0.97±0.14 


18.83 c 


0.006° 


6.3 e 


NGC 339 




-1.36±0.16 b 


-1.25±0.18 


18.75° 


0.04 e 


6.0 C 


Lindsay 38 




-1.59±0.10 b 


-1.52±0.12 


19.00 c 


0.04 e 


6.5 e 



a Iron abundance by Kraft & Ivans (2004). 

b Iron abundance in the ZW84 metallicity scale by Da Costa & Hatzidimitriou (1998, 
Kron 3, NGC 339) and Kayser ct al. (2007, Lindasy 38). 
c Iron abundance by C09. 

d Distance modulus and reddening by Ferraro et al. (1999, M79), Vargas et al. (2007, 
M55) and Dckany & Kovacs (2009, M53). The uncertainty on distance and reddening 
are 5% and 20%. 

"Distance modulus, reddening and age by Glatt ct al. (2009) using cluster isochroncs 
by Dotter et al. (2007). 
f Cluster ages: Salaris & Weiss (2002, M79, M55), Sarajedini (2009, M53). 
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Fig. 1. — Top - U-V,B-I Color-Color Diagram (C-CD) of Carina dSph. The curved box encloses candidate Carina stars (grey dots). 
The black dots represent probable field objects. Predicted distributions (Rocca-Volmerange et al. 2008) of elliptical (asterisks), spiral 
(pluses), SB (triangles), irregular and star burst (squares) galaxies as a function of the redshift are also showed. The vertical colored arrows 
mark the predicted (Castellani et al. 2002) peaks in B-V color for the halo (blue), the thick disk (green) and the thin disk (red). The 
black arrow shows the reddening vector. Bottom - Left - V, B-V Color-Magnitude Diagram (CMD) of Carina dSph. Bottom — Middle - 
V, B-V CMD of candidate field stars according to the C-CD selection. The error bars on the right represent the intrinsic photometric error 
in magnitude and in color. Rejection of thick-disk turnoff stars is obviously imperfect near £?-/~0.5, where the C-CD loses discriminating 
power. Bottom - Right — V,B-V CMD of candidate Carina stars according to the C-CD selection. A hint of the Carina HB near V«20.75 
and of the SGB at 22.5 < V < 23.75 confirm that discrimination is imperfect at these intermediate colors. The block of dots with B-I <0.5, 
V>23.5 is probably dominated by background galaxies rather than Carina turnoff stars. 
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M79 - NGC 1904 

[Fe/H] m =-1.64 



/n=15.61 
E(B-V)=0.01 



M55 - NGC 6809 

[Fe/H] KIM =-1.85 , 



fl= 13.69 
E(B-V)=0.09 




M53 - NGC 5024 

[Fe/H] lw =-2.02 



H= 16.24 
E(B-V)=0.01 




E(B-I)= 0.02 
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M | 
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Kron 3 

[Fe/HW=' 



-1.08 



ji= 18.83 
E(V-I)=0.008 

NGC ' 339 ' 

[Fe/HU=-1.36 




H= 18.75 
E(V-I)=0.06 




Lindsay 38 



[Fe/H] zl 



-1.59 



H= 19.00 
E(V-I)=0.05 
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Fig. 2.— Left column - the V, B-V CMD of the three old calibrating GCs, namely M79 (top), M55 (middle) and M53 (bottom). Middle 
column - same as the left, but for the V, B- 1 CMD. Right column - the V, V—I CMD of the three intermediate-age SMC calibrating 
clusters, namely Kron 3 (top), NGC 339 (middle) and Lindsay 38 (bottom). The red lines display the ridge lines, the red squares mark old 
HB stars and the red oval the contour of RC stars. The iron abundance, the distance modulus, and the reddening are also labeled. 
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-V [mog] B-l [mog] V-l [mog] 

Fig. 3. — Left column - comparison in the V, B-V CMD between candidate Carina stars (black dots) and the three old calibrating 
GCs, namely M79 (top), M55 (middle) and M53 (bottom). The distance modulus and the reddening adopted to overlap the SGB of old 
GCs with the SGB of Carina are labelled. Middle column - same as the left, but in the V, B-I CMD. Right column - comparison in the 
V, V-I CMD between Carina stars and the three intermediate-age SMC calibrating clusters, namely Kron 3 (top), NGC 339 (middle) and 
Lindsay 38 (bottom). The distance modulus and the reddening adopted to overlap the SGB of intermediate age SMC clusters with the 
SGB of Carina are labelled. 
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0.0 0.5 1.0 1.5 
B— V [mag] 

Fig. 4. — Top - cluster isochrones in the V, B-V CMD at fixed age (6 Gyr) and for a broad range of chemical compositions (BaSTI 
database). The vertical red dashed line shows the isochrone for an old (f3 Gyr), metal-poor ( [Fe/H] =-1 . 79 dex) stellar structure together 
with its ZAHB. The colored lozenges mark the RC color, while the red square marks the middle of the RR Lyrae instability strip 
(logT e =3.85). Bottom - cluster isochrones in the V, B-V CMD at fixed age (13 Gyr) and for a broad range of chemical compositions 
(BaSTI database). The vertical red dashed line shows the isochrone for an intermediate age (6 Gyr), metal-intermediate ([Fe/H]=-1.27 
dex) stellar structure together with its He-burning phases. The symbols are the same as in the top panel. 
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Fig. 5. — Left column - zoom across the RC and old HB region of Carina in the V, B-V (top), V, B-I (middle) and in the V, B-I 
(bottom) CMD. The red circle marks the position of the middle of the RR Lyrae instability strip and the red asterisk marks the position of 
the peak of RC stars. Right column - Predicted difference in B-V (top), in B-I (middle) and in V-I (bottom) colors between the middle 
of the RR Lyrae instability strip and the peak of RC stars as a function of iron abundance. Lines of different colors display predictions 
for different ages of intermediate-mass stars (see labeled values). The black asterisks mark the BaSTI grid points. The red dots display 
Carina, while the error bars the uncertainties affecting the difference in color and the metallicity estimate. 
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Old: IHUGyr [Fe/H] = -2.0±0. 1 dex 

Inter: 2.5H-4Gyr [Fe/H] = - 1 .5±0.1dex 
M = 20.15 
E(B-V) = 0.03 




Old: 11-M3Gyr [Fe/H] = -2.0±0.2dex 




Old: 11^-1 3Gyr [Fe/H] = -2.0±0.2dex 




0.0 0.5 1.0 1.5 

B-V [mag] 



Fig. 6. — Top - synthetic V, B-V CMD of Carina. The CMD was constructed by assuming two star formation episodes, with ages of 
11-^13 Gyr for the old and of 2.5-^4 Gyr for the intermediate-age stellar component (Monelli et al. 2003). The two star formation episodes 
include the same fraction of stars and have a mean metallicity of [Fe/H]=-1.5 and of -2.0 dex, respectively. The spread in metallicity is 
±0.1 dex for both the intcrmdiatc and the old sub-population. Intrinsic photometric error was added using a Gaussian with a mean equal 
to 0.02 mag. Middle - same as the top, but the synthetic CMD was constructed assuming a spread in metallicity of ±0.2 dex for the old 
and the intcrmdiatc-agc sub-population. Bottom — same as the middle, but the synthetic CMD was constructed assuming a spread in 
metallicity of ±0.5 dex for the intermdiate-age sub-population. 



